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C^) We propose a model which can be jointly calibrated to the corporate bond term 

structure and equity option volatility surface of the same company. Our purpose is to 
obtain explicit bond and equity option pricing formulas that can be calibrated to find 
a risk neutral model that matches a set of observed market prices. This risk neutral 
model can then be used to price more exotic, illiquid or over-the-counter derivatives. We 
observe that the model implied credit default swap (CDS) spread matches the market 
CDS spread and that our model produces a very desirable CDS spread term structure. 
This is observation is worth noticing since without calibrating any parameter to the CDS 
spread data, it is matched by the CDS spread that our model generates using the available 
I information from the equity options and corporate bond markets. We also observe that 

our model matches the equity option implied volatility surface well since we properly 
account for the default risk premium in the implied volatility surface. We demonstrate 
the importance of accounting for the default risk and stochastic interest rate in equity 



u 
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£ — . option pricing by comparing our results to Fouque et al. (20031, which only accounts for 

stochastic volatility. 

Keywords: Credit Default Swap, Defaultable Bond, Defaultable Stock, Equity Options, 
Stochastic Interest Rate, Implied Volatility, Multiscale Perturbation Method. 



Contents 



1 Introduction [2] 

'Department of Mathematics, University of Michigan, Ann Arbor, MI 48109, USA, e- 
mail: {erhan, boy} @umich.edu. 

^We would like to thank the two anonymous referees and the anonymous AE for their constructive com- 
ments, which helped us improve our paper. 

*This work is supported in part by the National Science Foundation, under grant DMS-060449. 



1 



2 A Framework for Pricing Equity and Credit Derivatives [5] 

2.1 The model 

2.2 Equity and credit derivatives [7] 

3 Explicit Pricing Formulas for Credit and Equity Derivatives [9] 

3.1 Pricing equation 

3.2 Asymptotic expansion [10] 

3.3 Explicit pricing formula [T3] 

4 Calibration of the Model [TH] 

4.1 Data description [19] 

4.2 The parameter estimation [2U] 

4.3 Model implied CDS spread matches the observed CDS spread [22] 

4.4 Fitting Ford's implied volatility [22] 

4.5 Fitting the implied volatility of the index options [23] 



1 Introduction 



Our purpose is to build an intensity-based modeling framework that can be used in trading 
and calibrating across the credit and equity markets. The same company has stocks, stock 
options, bonds, credit default swaps on these bonds, and several other derivatives. When this 
company defaults, the payoffs of all of these instruments are affected; therefore, their prices 
all contain information about the default risk of the company. 

We build a model that can be jointly calibrated to corporate bond prices and stock 
options, and can be used to price more exotic derivatives. In our framework we use the 
Vasicek model for the interest rate, and use doubly stochastic Poisson process to model the 
default of a given company. We assume that the bonds have recovery of market value and 
that stocks become valueless at the time of default. Using the multi-scale modeling approach 



of Fouque et al. (20031 we obtain explicit bond pricing equation with three free parameters 
which we calibrate to the corporate bond term structure. On the other hand, stock option 
pricing formula contain seven parameters, three of which are common with the bond option 
pricing formula. (The common parameters are multiplied with the loss rate in the bond 
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pricing formula.) We calibrate the remaining set of parameters to the stock option prices. 
This hybrid model, therefore, is able to account for the default risk premium in the implied 
volatility surface. 

The calibration results reveal that our model is able to produce implied volatility surfaces 
that match the data closely. We compare the implied volatility surfaces that our model 
produces to those of Fouque et al. ( 2003[ ). We see that even for longer maturities our model 
has a prominent skew: compare Figures [5] and [6] Even when we ignore the stochastic volatility 
effects, our model fits the implied volatility of the Ford Motor Company well and performs 
better than the model of Fouque et al. (2003); see Figure [4j This points to the importance 
of accounting for the default risk for companies with low ratings. 

Once the model parameters are calibrated, the model can be used to compute the prices 
of more exotic options. To test whether our model produces correct prices we use the CDS 
spread data and show that the model implied CDS spread matches the "out of sample" CDS 
data. To compute the CDS spread, under our assumption on the recovery, one needs to 
reconstruct the term structure of the treasury and the corporate bonds. Moreover, one needs 



to separate the loss rate from the other parameters in the bond pricing formula (see (2.14) or 



(4.7) for the CDS spread formula). This separation is possible since we calibrate our model 
to corporate bond data and stock option data jointly as described above. The model-implied 
CDS spread time series matches the observed CDS spread time series of Ford Motor Company 
for over a long period of time; see Figures [T] and [2j This is an interesting observation since 
we did not make use of the CDS spread data in our calibration. This observation also shows 
that one can use our model to trade across different markets that contain information about 
the default risk of a given firm. 

Our model has three building blocks: (1) We model the default event using the multi-scale 



stochastic intensity model of Papageorgiou and Sircar (2008). We also model the interest rate 



using an Ornstein-Uhlenbeck process (Vasicek model). As it was demonstrated in Papageor- 



giou and Sircar (2008), these modeling assumptions are effective in capturing the corporate 



yield curve; (2) We take the stock price process to follow a stochastic volatility model which 
jumps to zero when the company defaults. This stock price model was considered in iBa yrak^l 



tar 



(2008). Our model specification for the stock price differs from the jump to default models 



for the stock price considered by Carr and Linetsky (2006) and Linetsky (2006), which take 



the volatility and the default intensity to be functions of the stock price; (3) We also account 
for the stochastic volatility in the modeling of the stocks since even the index options (when 
there is no risk of default) possess implied volatility skew. We model the volatility using 



the fast scale stochastic volatility model of Fouque et al. (2000). We demonstrate on index 



options (when there is no risk of default) that (see Section 4.5), we match the performance 
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of the two time scale volatility model Fouque et al. (20031. The latter model extends Fouque 



et al. ( 2000 1 by including a slow factor in the volatility to get a better fit to longer maturity 



option. We see from Section 4.5 that when one assumes the interest rate to be stochastic, the 



calibration performance of the stochastic volatility model with only the fast factor is as good 
as the two scale stochastic volatility model, which is why we choose the volatility to be driven 
by only the fast factor. Even though the interest rate is stochastic in our model, we are able 
to obtain explicit asymptotic pricing formulas for stock options. Thanks to these explicit 
pricing formulas the inverse problem that we face in calibrating to the corporate bond and 
stock data can be solved with considerable ease. Our modeling framework can be thought of 



as a hybrid of the models of Fouque et al. (2000), which only considers pricing options in a 



stochastic volatility model with constant interest rate, and Papageorgiou and Sircar (2008), 



which only considers a framework for pricing derivatives on bonds. Neither of these models 
has the means to transfer information from the equity markets to bond market or vice versa, 
which we are set to do in this paper. We should also note that our model also takes input 
from the treasury yield curve, historical stock prices, and historical spot rate data to estimate 
some of its parameters (see Section El) . 



Our model extends Bayraktar (2008) by taking the interest rate process to be stochastic, 



which leads to a richer theory and more calibration parameters, and therefore, better fit 
to data: (i) When the interest rate is deterministic the corporate bond pricing formula 
turns out to be very crude and does not fit the bond term structure well (compare (2.57) 
in Bayraktar (2008) and (4.1)); (ii) With deterministic interest rates the bond pricing and 



the stock option pricing formulas share only one common term, "the average intensity of 
default" (this parameter is multiplied by the loss rate in the bond pricing equation, under 
our loss assumptions). Therefore, the default premium in the implied volatility surface is not 
accounted for as much as it should be. And our calibration analysis demonstrates that this 
has a significant impact. When the volatility is taken to be constant, both our new model 



and the model in Bayraktar (2008) have three free parameters. The model in Bayraktar 



(2008) produces a below par fit to the implied volatility surface (see e.g. Figure 5 in that 
paper), whereas our model produces an excellent fit (see Section [4~4| and Figure [4]); (iii) To 
calculate the CDS spread, in the constant interest rate model, one needs to separate the loss 
rate and the average intensity of default. This is again established calibrating the model to 
the bond term structure data and the stock option implied volatility surface. The estimates 



for the average intensity and the loss rate are not as accurate in Bayraktar (2008) as it is in 



our model because of (i) and (ii). This crude estimation leads to a poor out of sample match 
to the CDS spread time series. 



The other defaultable stock models are those of Carr and Linetsky (2006), Linetsky (2006 ) 



and Carr and Wu (2006), which assume that the interest rate is deterministic. Carr and 
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Linetsky ( 2006 1 , Linetsky ( 2006 1 take the volatility and the intensity to be functions of the 



stock price and obtain a one-dimensional diffusion for the pre-default stock price evolution. 
Using the fact that the resolvents of particular Markov processes can be computed explicitly, 



they obtain pricing formulas for stock option prices. On the other hand Carr and Wu (20061 



uses a CIR stochastic volatility model and also models the intensity to be a function of 
the volatility and another endogenous CIR factor. The option prices in this framework are 
computed numerically using inverse the Fourier transform. We, on the other hand, use 
asymptotic expansions to provide explicit pricing formulas for stock options in a framework 
that combines a) the Vasicek interest rate model, b) fast-mean reverting stochastic volatility 
model, c) defaultable stock price model, d) multi-scale stochastic intensity model. 



Our calibration exercise differs from that of Carr and Wu (2006) since they perform a 



time series analysis to obtain the parameters of the underlying factors (from the the stock 
option prices and credit default swap spread time series), whereas we calibrate our pricing 
parameters to the daily implied volatility surface and bond term structure data. Our purpose 
is to find a risk neutral model that matches a set of observed market prices. This risk neutral 
model can then be used to price more exotic, illiquid or over-the-counter derivatives. For 



further discussion of this calibration methodology we refer to Cont and Tankov (2004) (see 



Chapter 13), Fouque et al. (2000), Fouque et al. (2003) and Papageorgiou and Sircar (2008). 
We also provide daily prediction of the CDS spread only using the data from the bond term 
structure and implied volatility surface of the options. 

The rest of the paper is organized as follows: In Section 2, we introduce our modeling 
framework and describe the credit and equity derivatives we will consider and obtain an 
expression for the CDS spread under the assumption that the recovery rate of a bond that 
defaults is a constant fraction of its predefault value. In Section 3, we introduce the asymp- 
totic expansion method. We obtain explicit (asymptotic) prices for bonds and equity options 
in Section 3.3. In Section 4, we describe the calibration of our parameters and discuss our 
empirical results. Figures, which show our calibration results, are located after the references. 



2 A Framework for Pricing Equity and Credit Derivatives 
2.1 The model 

Let (fl,Tt,F) be a complete probability space supporting (i) correlated standard Brownian 
motions W t = (W?, W t \ W? , Wf, W t 4 ), t > 0, with 



E[W?, Wt] = Pi t, E[Wi, Wi] = Pij t, i,j e {1, 2, 3, 4}, t > 0, 



(2.1) 
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for some constants pi,pij G (—1, 1), and (ii) a Poisson process N independent of W . Let us 
introduce the Cox process (time-changed Poisson process) N t = N(f* \ s ds), t>0, where 

X t = f(Y t ,Z t ), 

dYt = \m-Y t )dt+ U -^dWl Y = y, (2.2) 
e ye 

dZ t = Sc(Z t )dt + V6g(Z t )dW?, Z = z, 

in which e, 5 are (small) positive constants and / is a strictly positive, bounded, smooth 
function. We also assume that the functions c and g satisfy Lipschitz continuity and growth 
conditions so that the diffusion process for Z t has a unique strong solution. We model the 
time of default as 

t = inf{t > : N t = 1}- (2.3) 
We also take interest rate to be stochastic and model it as an Ornstein-Uhlenbeck process 

dr t = (a- (3r t )dt + r]dW}, r = r, (2.4) 

for positive constants a, (3, and r\. 

We model the stock price as the solution of the stochastic differential equation 

dX t = X t (r t dt + atdWt ~ d {^ t ~ Xudu ) ) ' ^ = x, (2.5) 
where the volatility is stochastic and is defined through 

a t = a(Y t ); dY t = (^(m - Y t ) - V -^K{Y t )^ dt + ^0 = V- (2.6) 

Here, A is a smooth, bounded function of one variable which represents the market price of 
volatility risk. The function a is also a bounded, smooth function. Note that the discounted 
stock price is a martingale under the measure P, and at the time of default, the stock price 
jumps down to zero. The pre-banktruptcy stock price coincides with the solution of 

dX t = {r t + \ t )X t dt + a t X t dW t °, X = x. (2.7) 



It will be useful to keep track of different flows of information. Let F = {J~t,t > 0} be the 
natural filtration of W. Denote the default indicator process by It = l{ T <t}) t > 0, and let 
I = {2 t , t > 0} be the filtration generated by I. Finally, let G = {Gt, t > 0} be an enlargement 
of F such that Q t = T t \l It, t > 0. 
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Since we will take e and 5 to be small positive constants, the processes Y and Y are 



fast mean reverting, and Z evolves on a slower time scale. See Fouque et al. (2003) for 



an exposition and motivation of multi-scale modeling in the context of stochastic volatility 
models. 



We note that our specification of the intensity of default coincides with that of Papa- 



georgiou and Sircar (2008), who considered only a framework for pricing credit derivatives. 



Our stock price specification is similar to that of Linetsky (2006) and Carr and Linetsky 



(2006) who considered a framework for only pricing equity options on defaultable stocks. 



Our volatility specification, on the other hand, is in the spirit of Fouque et al. (2000) 



Bayraktar (2008) considered a similar modeling framework to the one considered here, 



but the interest rate was taken to be deterministic. In this paper, by extending this modeling 
framework to incorporate stochastic interest rates, we are able to consistently price credit 
and equity derivatives and produce more realistic yield curve and implied volatility surfaces. 
We are also able to take the equity option surface and the yield curve data as given and 
predict the credit default swap spread on a given day. Testing our model prediction against 
real data demonstrates the power of our pricing framework. 



2.2 Equity and credit derivatives 



In our framework, we will price European options, bonds, and credit default swaps of the 
same company in a consistent way. 

1. The price of a European call option with maturity T and strike price K is given by 

i-T 



C(t; T,K)=M 



exp 



r s ds (X T - K) + l {T>T} 



exp f- J (r s + \ s )ds} (X T - K) 



Ft 



(2.8) 



in which the equality follows from Lemma 5.1.2 of Bielecki and Rutkowski (2002). (This 



lemma, which lets us write a conditional expectation with respect to Qt in terms of conditional 
expectations with respect to Tt, will be used in developing several identities below). Also, 



see Linetsky (2006) and Carr and Linetsky (2006) for a similar computation 
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On the other hand, the price of a put option with the same maturity and strike price is 



Put (t;T) = E 



expl-J^ r s ds j (K - X t ) + 1 {t>t} Gt + E exp ( - J r s ds j K 1{ T <T} \G; 



= l{r> t} E 



exp 



(r s + X s )ds {K-X T y 



KE 



exp 



r s ds 



T 



KE 



exp 



T 



(r s + X s )ds 



T 



(2.9) 



2. Consider a defaultable bond with maturity T and par value of 1 dollar. We assume the 
recovery of the market value, introduced by Duffie and Singleton (19991. In this model, if 
the issuer company defaults prior to maturity, the holder of the bond recovers a constant 
fraction 1 — I of the pre-default value, with / 6 [0, 1]. The price of such a bond is 



B c {t;T) = E 



E 



exp 



r s ds l{ r>T } +exp 



ds l {r < r} (l-0# c (T-;T) 



exp f- jf (r s + l X s )ds^j 



T, 



(2.10) 



on {r > t}, see Duffie and Singleton ( |1999 ) and Schonbucher (1998 1 . 



3. Consider a credit default swap (CDS) written on B c , which is a insurance against 
losses incurred upon default from holding a corporate bond. The protection buyer pays a 
fixed premium, the so-called CDS spread, to the protection seller. The premium is paid on 
fixed dates T = (T 1; • • • , Tm), with Tm being the maturity of the CDS contract. We denote 
the CDS spread at time t by c ds (t;T). Our purpose is to determine a fair value for the CDS 
spread so that what the protection buyer buyer is expected to pay, the value of the premium 
leg of the contract, is equal to what the protection seller is expected to pay, the value of 
the protection leg of the contract. For a more detailed description of the CDS contract, see 



Bielecki and Rutkowski 


(2002 


1 or 


Schonbucher 


(2003) 



The present value of the premium leg of the contract is 



M 

^exp 



m=l 



Premium^; T) = c ds (t; T) E 

= l {T>t} c ds (t;T)J2^ 



M 



m=l 



exp 



r s ds l {T>rm} 



(r s + X s )ds 



(2-11) 



T 



8 



in which we assumed that t < T±. The present value of the protection leg of the contract 
under our assumption of recovery of market value is 



Protection^; T) = l{ r>t }E 



eXP (~/ r sds) 1{ t <t m }IB c (t-;T m ) 



(2.12) 



Adding (2.10) and (2.12), we obtain 



Protection^; T) + B c (t; T M ) = E 



exp 



r s ds l{ r >r} +exp 



l{r>t}E 



exp [ — j r s ds 



Ft 



ds\\ {T < T} B c {r--T) 



(2.13) 



where the last equality is obtained by setting I = in (2.10). 



Now, the CDS spread can be determined, by setting Protection(t; T) = Premium(i; T) 



and using equations (2.11) and (2.13), as 



E 



c ds (t;T) = l {T>t} - 



exp ( - Jt M r s ds 



F t 



E 



exp(-J^ M (r s + l\ s )ds 



M 



m=l 



exp 



(r s + X s )ds 



F t 



(2.14) 



3 Explicit Pricing Formulas for Credit and Equity Derivatives 



3.1 Pricing equation 

Let P e ' 5 denote 

P e > s (t,X t ,r t ,Y t ,Y t ,Z t ) = E 



exp 



(r s + l\ s )ds) h(X T ) 



F, 



(3-1) 



When / = 1 and h(X T ) = (X T -K)+, P^ 5 is the price of a call option (on a defaultable 
stock). On the other hand, when h{Xx) = 1, P e ' S becomes the price of a defaultable bond. 



Using the Feynman-Kac formula, we can characterize P t,s as the solution of 

£ e > s P e > 5 (t,x,r,y,y,z) = 0, 
P e ' 5 (T,x,r,y,y,z) = h(x), 



(3.2) 
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where the partial differential operator C e ' is defined as 



e Je V e 



(3.3) 



in which 



A 2 

dy 2 



d 2 



8 2 r- d 2 ^ d 2 _ ,- d 2 ^8 



dydy ' 



A = P2(r{y)vy/2x—-^- + pi 2 i]^V2—— + P4,cr(y)i>V2x--— + PuW^^-^ ~ Mv)"^^- 
oxay oroy oxay oroy ay 



did 2 d d d 2 1 d 2 



Mi 4 aWwWx^- + VPi39(z)^g- z , M 2 4 + ^ 2 (z) J^, 

-M 3 - /023^\/25(z)^-^- + ,0341^/25(2)^^-. 



3.2 Asymptotic expansion 

We construct an asymptotic expansion for P e,s as e, 5 — > 0. First, we consider an expansion 
of P t ' in powers of V5 

P^ s = P' + y/spf + SP* + ... (3.4) 



By inserting (3.4) into (3.2) and comparing the 5° and 5 terms, we obtain that Pq satisfies 

1 



^c + l=c 1 + c 2 ) rr, -0. 



Po( T iX,r,y,y,z) = h(x), 



(3.5) 



and that P{ satisfies 



1 1 -*-2 I -»l 

P[{T,x,y,y,z,r) = 0. 



Mt + ^=M 3 ) Ply 



(3.6) 



Next, we expand the solutions of (3.5) and (3.6) in powers of y/e 



P e Q = P + + eP 2 , + e 3 / 2 P3,o + • • • 

Pi = Po,i + V~tPi,i + eP 2 ,i + e 372 ^,! + • • 



(3.7) 
(3.8) 
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Inserting the expansion for Pq into (3.5) and matching the 1/e terms gives CqPq = 0. We 
choose Pq not to depend on y and y because the other solutions have exponential growth at 



infinity (see e.g. Fouque et al. (2003)). Similarly, by matching the l/y/e terms in (3.5) we 



obtain that -Co-Pi.o + £iPo = 0. Since C\ takes derivatives only with respect to y and y, we 
observe that £o-Pi,o = 0. We choose P^o not to depend on y and y. 



Now equating the order-one terms in the expansion of (3.5) and using the fact that 
£i-Pi,o = 0, we get that 

£ P 2 ,o + £2Po = 0, (3.9) 



which is a Poisson equation for P 2 fi (see e.g. Fouque et al. (2000)). The solvability condition 
for this equation requires that 

(C 2 )P = 0, (3.10) 

where (•) denotes the averaging with respect to the invariant distribution of (Y t ,Y t ), whose 
density is given by 



l 



l 



2-kvv 
Let us denote 



- exp 



2(1 -p| 4 ) 



y — m 



v 



+ 



y — m 



v 



2^24 



(y - m)(y - fh) 



To demonstrate the effect of averaging on £2, let us write 



(A) : =| + ^ 2 V^ + (r + A(z))x^ + (a-/?r)|- + a lW ^- r , , /; ^ 



+ n 7 ? 



1 2 <9 2 



Together with the terminal condition 



P (T,x,r,z) = h(x), 



(3.11) 
(3.12) 



(r + l\(z))- 
(3.13) 



(3.14) 



equation (|3.10|) defines the leading order term Pq. On the other hand from (3.9), we can also 

^2,0 = -C \C 2 - (C 2 ))P . (3.15) 



deduce that 



Matching the y/e order terms in the expansion of (3.5) yields 



£ -P 3 ,o + £1^2,0 + £ 2 Pi, = 0, 



(3.16) 



which is a Poisson equation for i^o- The solvability condition for this equation requires that 
(£ 2 P li0 ) = -<AP 2 ,o> = (CiC \C 2 - <£ 2 »>P , (3.17) 
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which along with the terminal condition 



Pi fi (T,x,r,z) = 0, 



(3.18) 



completely identifies the function Pi o- To obtain the second equality in (3.17) we used (3.15 ). 



Next, we will express the right-hand side of (3.17) more explicitly. To this end, let tp, k, 
and <p be the solutions of the Poisson equations 



AWKi/) = cr(y) - 0"! ^o^(y) = <? 2 {y) - o{, and C (j)(y, z) = (f(y, z) - \(z)), (3.19) 



respectively. First observe that 

1 B 2 P B 2 P / BP 

(C 2 -{C 2 ))P = 2( ff 2 te)-^V 2 ^+Wy)-^)w^+i(/(^)-A(^)) [x^ 



-Po 
(3.20) 



Now, along with (3.19), we can write 

12 



c \c 2 - (c 2 ))p = \<y)* 2d ^!- + mm*^ + 1 ttv, *) - ft) • (3.21) 

Applying the differential operator £1 to the last expression yields 

,— BP — B f BP 

{C x Cu X {C 2 - (C 2 )))P = lp 2 vV2{cr ( t }y )x 2 -^ + lp 12W V2(cj )y )— ix-^ - P 



1 



V2 -(any 



d ( 2 d 2 P 



x-— x 



dx \ Bx 2 



B ( B 2 P 



+ (aipy)rjpix— x 



Bx \ BxBr 



+ pi 4 r]i)V2 ( -(ny) 



B ( 3 2 P ( 



1. ; ( r 2')' 2 l\ 

2 Or v <9x 2 / ^ /ri 9r V~ BxBr J 



-^(liA^x^ + iA^ 



Bx 2 



3 2 ft\ 



•)y)rjpix 



BxBr I 



(3.22) 



Finally, we insert the expression for Pf in (3.8) into (3.6) and collect the terms with the 



same powers of e. Arguing as before, we obtain that Po 5 i is independent of y and y and 
satisfies: 



(£ 2 )P ,i = -(Mi)P , P ,i(T,x) = 0. 



(3.23) 
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3.3 Explicit pricing formula 



We approximate P e ' defined in (3.1) by 



P e ' s = P + ViPi,o + VSP i . 



(3.24) 



Since the Vasicek interest rate process is unbounded, which implies that the potential term 



in £2 01 the discounting term in (3.1) is unbounded, the arguments of Fouque et al. (2003) 



can not be directly used. However as in Cotton et al. (2004) and Papageorgiou and Sircar 



(2008), one can write 



P e >\t, X t , n, Y t , Y t , Z t ) = B(t, T)E 



in which 



T 



exp 



T 



lX s ds) h(X T ) 



T exp ( — f Q r s ds 



dF 



s(o, r) 



and 



B(t,T) =E 



exp 



T 



r*ds 



B(t,T)F*> s (t,X t ,r t ,Y t ,Y t: 
(3.25) 

(3.26) 
(3.27) 



Now, the analysis of Fouque et al. (2003) can be used to approximate F e,s (t, x,r,y,y, z). As a 
result of this analysis for each (t, x, r, y, y, z), there exists a constant C such that \P e ' S — P e ' S \ < 
C ■ (e + 5) when h is smooth, and \P £ ' S — P e,5 | < C ■ (elog(e) + S + \fe$) when h is a put or a 
call pay-off. In what follows, we will obtain Pq, P^o and Pq 1 explicitly. 



Our first objective is to develop a closed-form expression for Po> the solution of (3.10) 



and (3.14). 



Proposition 3.1. The leading order term Po in (3.24) ^ s given by: 

/oo 
h(exp(u)) 
-00 

where 



1 



, (u - m(t,T)) 2 . , 
exp(- V )du, (3.28) 



2v(t,T) 



B$(t, r; z, T, I) = exp ( - l\(z)(T - t) + a(T - t) - b(T - t)r), 
in which the functions a(s) and b(s) are defined as: 



a(s) 



(3.29) 



(3.30) 
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and b(s) = (1 — exp(— [3s))/ (3. On the other hand, 



2 2r ] p 1 a 1 
v t,T = I cr 2 + 



P P 2 



2r]pi(Ji 2rf_ 

P 2 pz 



exp(-P(T-t)) 



V „, (2r\p\o~\ 3rf 



2P 



3 exp(-2P(T-t)) 



P 2 + 2/?V 



and 



mt, T = log(x) + A • (T - t) - a(T - 1) + b(T - t)r - -v{t, T). 



Proof. By applying the Feynman-Kac theorem to (3.10) and (3.14) we have that 

rT 



Po(t, x, z, r) = E exp 
where the dynamics of S is given by 



(r s + fX{z))ds^j h{S T ) 



S t = x,r t = r 



dS t = (r t + \{z))S t dt + a 2 S t dW?, 

in which W° is a Wiener process whose correlation with W 1 is p\ = f^pi- 
Let us define 



P (t,x,r)=E 



exp 



T \ ~ 
r s ds \ h(S T ) 



S t = x,r t = r 



in which 



Then 



dS t = r t S t dt + a 2 S t dW? . 



P (t, x, z, r) = e- ,X « ( T - f )Po(t,^exp(A( 2 )(T - t)), r). 



forward measure P T through the Radon-Nikodym derivative (3.26) 



in which 



F t 



B(t,TY 



(3.31) 



(3.32) 



(3.33) 



(3.34) 



(3.35) 

(3.36) 
(3.37) 



Now, by following Geman et al. (1995), we change the probability measure P to the 



We can obtain the following representation of Pq using the T-forward measure 

P (t,S t ,r t ) = B(t,T)E T \h(S T )\Tt] = B(t, T)E T [h{F T )\T t ] , (3.38) 



(3.39) 



which is a P T martingale. Note that an explicit expression for B(t, T) is available since rt is 
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a Vasicek model, and it is given in terms of the functions a and b: 



B(t, T) = exp(a(T - t) - b(T - t)r t ). 



By applying Ito's formula to (3.39), we observe that the dynamics of F are 



dF t = FtiaxdW? + b(T - t)ridW}) 



(3.40) 



(3.41) 



in which W l is a F T Brownian motion whose correlation with the W° (which is still a 
Brownian motion under F T ) is p\. Given Xt and B(t,T), the random variable log-F^ is 
normally distributed with variance 



Vt T = a\{T -t)+r] 2 / b 2 (T - s)ds + 2r]p x a 2 \ b(T - s)ds 
Jt Jt 



2 2r 1Pl a 2 ,rf UT ,\ , 



(3 2 



exp(-/3(T-i)) 



(3.42) 



and mean 



m(t,T)=logF t -^£ (a 2 2 + b 2 (T-s)r ] 2 + p 1 a 2 b(T-s)r ] )ds = log ^ 

Now the result immediately follows. 

An immediate corollary of the last proposition is the following: 



1 



£>(i, T) / 2 



t*,T. (3.43) 



□ 



Corollary 3.1. i) When I = 1, h(x) = (x — K) + , then (3.28) becomes 



C (t, x, z, r) = xN{d x ) - KB c (t, r; z, T, l)N(d 2 ), 
in which N is the standard normal cumulative distribution function and 

l °£ KB5(t*r-,z,T,l) ± h V t,T 



(3.44) 



/VtT 



(3.45) 



ii) When I = 1, and h{x) = {K — x) + , then (3.28) becomes 



Puto(t, x, z, r) = -x + xN(di) - KBffrt, r; z, T, l)JV(cfc) + KB%(t, r; z, T, 0). (3.46) 



Hi) When h(x) = 1, then (3.28) coincides with (3.30) in Papageorgiou and Sircar (2008) 
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Proposition 3.2. The correction term \feP\Q is given by 



d 2 P dP 



dxda da 



dx 2 d- 



a 



d 2 P 

drjdx 



d 2 P 
dxda ' 



(3.47) 



in which 

V{ = ^/e (l p 2 i/V2(a(f> y ) - uV2^(AKy)^j , V 2 = \^p A v\Fi(pKy) , 

V£ = Ve{pi2W^(4>y)), K = ~Ve (^piAr]i)V2(Ky) - p A vV2{a^ y )rip 1 + piAwV2(ipy)a l pf S j , 

V$ = -^/e(p u r]i>V2(ipy) px) , Vq = \Ze(-pii>V2{aip^)r)pi + pwqDV^^^axp 2 - vV2(Aif}y)r]pi). 

(3.48) 



Proof. Recall that P\^ is the solution of (3.17) and (3.18) and that the right-hand-side of 



(3.17) is given by (3.22). The result is a simple algebraic exercise given the following four 



observations: 

1) x n -j^ commutes with (£2)- 

2) -(T-t)(x n £,)P solves: 



(C 2 )u = x 



dx % 



P , u(T,x,r; z) = 0. 



(3.49) 



3) By differentiating (3.14) with respect to a, we see that — jf®- also solves 



{C 2 )u 



m 

dr 



u(T, x, r; z) = 0. 



(3.50) 



4) Using 1) and 2) above and the equation we obtain by differentiating (3.10) with respect 
to 77, we can show that 1 / 77 • {°^iPi x Q x ^° a ~ 7^) solves 



(C 2 )u 



d 2 Pp 
dr 2 



u(T, x, r; z) = 0. 



(3.51) 



□ 



Remark 3.1. By differentiating (3.10) with respect to r, we obtain 

dP 



(A 



dr 



ox dr 



(3.52) 
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Using observation 2 in the proof of Proposition 
solves 



3.2 



we see that 4 



i (-(T-t)(x 



dPg 

dx 



{C 2 )u 



m 

dr 



, u(T, x, r; z) = 0. 



Now, it follows from observation 3 in the proof of Proposition \3.2\ that 

dP 1 



da (3 



rr i \ I dP ° p\^ 9P °\ 



Using this identity, we can express ( 3.47) only in terms of the "Greeks" 



Pc) + 



dPp 
dr 



(3.53) 



(3.54) 



Next, we obtain an explicit expression for Pq,i> the solution of (3.23). We need some 



preparation first. By differentiating (3.10) with respect to z, we see that solves 



- OP - 
{Ct)u = -X'(z)x-^ + / X'(z)P , u(T, x, r; z) = 0. 



As a result (see Observation 2 in the proof of Propostion 3.2) 

dP ° {T -t)X'(z)(x d ^-lP 



dz \ dx 

from which it follows that — {M\)Pq can be represented as 



Proposition 3.3. The correction term V6Pq i is given by 



(T-t) 



2 ' 2 d 2 p 



dx 2 



" Ix^ + lPo) -(T-t)(x 



Dx 



d 2 P , dP 
da 
dP 



dadx 



d 2 p 







drdx d 



in which 



Vf = VsX'(z)a 1 p 3 g(z), V* = VSX '(z) V p 13 g(z). 



(3.55) 



(3.56) 



dx J \ dxdr dr 

(3.57) 



(3.58) 



(3.59) 



Proof. We construct the solution from the following observations and superposition since 
(£2) is linear: 



1) We first observe that ^ T (x n Jr-jt)Po solves 



0" 



(C 2 )u = — (T - t) ( x n — ) P Q , u(T, x, r;z) = 



(3.60) 
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2) Next, we apply (£2) on (T — t)^- and obtain 



<£ 2 ) (T - 1) 



dr 



— — + {T-t)[ -x^— + (3—- + P 
or V ox or 



as a result of which we see that 

dP {T-t) 2 ( dP { 



1 

P 



solves 



da 



(jC 2 )u = {T-t) 



x— Pq) +{T-t) 



dx 



dr 



m 

dr 



u(T, x, r; z) = 0. 



(3.61) 



(3.62) 

(3.63) 

□ 



4 Calibration of the Model 



In this section, we will calibrate the loss rate I and the parameters 



{A, V{, Vi, VI VI, Vi, Vi, Vf, V 2 5 }, 



which appear in the expressions (3.28), (3.47), and (3.58) on a daily basis (see, e.g., Fouque 



et al. ( 2003 ) and Papageorgiou and Sircar ( 2008 ) for similar calibration exercises carried out 



only for the option data or only for the bond data). We demonstrate this calibration on Ford 
Motor Company. Note that there are some common parameters between equity options and 
corporate bonds. Therefore, our model will be calibrated simultaneously to both of these 
data sets. We will also calibrate the parameters of the interest rate and stock models to the 
yield curve data, historical spot rate data and historical stock price data. Next, we test our 
model by using the estimated parameters to construct an out-of-sample CDS spread time 
series (3 year and 5 year), which matches real quoted CDS spread data over the time period 
(1/6/2006 - 6/8/2007) quite well. 

We also look at how our model-implied volatility matches the real option implied volatility. 



We compare our results against those of Fouque et al. (2003). We see that even when we 



make the unrealistic assumption of constant volatility, our model is able to produce a very 
good fit. 

Finally, in the context of index options (when A = 0), using SPX 500 index options data, 
we show the importance of accounting for stochastic interest rates by comparing our model 



to that of Fouque et al. (2000 2003). 
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4.1 Data description 

• The daily closing stock price data is obtained from nnance.yahoo.com. 

• The stock option data is from OptionMetrics under WRDS database, which is the same 



database used in Carr and Wu (2006) 



— For index options, SPX 500 in our case, we use the data from their Volatility 
Surface file. The file contains information on standardized options, both calls and 
puts, with expirations of 30, 60, 91, 122, 152, 182, 273, 365, 547, and 730 calender 
days. Implied volatilities there are interpolated data using a methodology based 
on kernel smoothing algorithm. The interpolated implied volatilities are very close 
to real data because there are a great number of options each day for SPX 500 
with different maturities and strikes. The calibration results for index options are 
presented in Figure [7] and only the data set on the June 8, 2007 is used. 

— On September 15, 2006 (Friday) Ford announced that it would not be paying divi- 
dends (see e.g. , http:/ /money. cnn.com / 2006 /09/15/ news / companies / ford/index. htm) 
Therefore, call options on Ford after that date do not have early exercise premium 
starting from Sep 18, 2006. We use Ford's implied volatility surface data from 
9/18/2006 to 6/8/2007. when we consider Ford Motor Company's options. We 
excluded the observations with zero trading volume or with maturity less than 9 
days. The calibration results are used to construct Figures [T]|6j In particular, the 
implied volatility surface data from 9/18/2006 to 6/8/2007 is used to construct [2] 

As opposed to the options on the index there are not as many individual company 
options; and we find that the results given by using interpolated implied volatilities 
in the Volatility Surface File and data implied volatilities differ. This may be due 
to the fact that there are a limited number of option prices available for individual 
companies; i.e., there may not be enough data points for the implied volatilities to 
be accurately interpolated. Therefore, we use the Option Price file, which contains 
the historical option price information, of the OptionMetrics database 

For each day we U.S government Treasury yield data with maturities: 1 month, 3 monts, 
6 months, 1 year, 2 years, 3 years, 5 years, 7 years, 10 years, 20 years. This data set 
is available at: www.treasury.gov/offices/domestic-finance/debt-management/interest- 
r ate / yield . shtml . 

Corporate bond and CDS data is obtained from Bloomberg. Number of available bond 
quotes and bond maturities vary. Typically there are around 15 data points, for exam- 
ple, on June 8th, we have the following maturities: 0.60278,1.0222,1.1861,1.3139,1.4083, 
1.5944, 2.3889, 2.6028, 3.0194, 3.2694, 3.3972, 3.6472, 4.1722, 4.3806, 6.3139, 9.5194. 
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4.2 The parameter estimation 



The following parameters can be directly estimated from the spot-rate and stock price his- 
torical data: 



1. The parameters of the interest rate model {a, f3, 77} are obtained by a least-square fitting 



to the Treasury yield curve as in Papageorgiou and Sircar (2008). 



2. p\ = f^pi, the "effective" correlation between risk- free spot rate r (we use the one- 



month treasury bonds as a proxy for r) and stock price in (3.34) is estimated from 
historical risk-free spot rate and stock price data. 



3. o"2, the "effective" stock price volatility in (3.34) is estimated from the historical stock 
price data. 

Now, we detail the calibration method for I, X(z) and {V{, Kf, Kf, V^, V^, V£, Vf, V$}. We 
will minimize the in-sample quadratic pricing error using non-linear least squares to calibrate 
these parameters on a daily basis. This way we find a risk neutral model that matches a set 
of observed market prices. This risk neutral model can then be used to price more exotic, 
illiquid or over-the-counter derivatives. This practice is commonly employed; and for further 



discussion of this calibration methodology we refer to Cont and Tankov (2004) (see Chapter 
13 and the references therein). 

Our calibration is carried out in two steps in tandem: 

Step 1. Estimation of IX and {IV%, IV2} from the corporate bond price data. 

The approximate price formula in ( 3.24[ ) for a defaultable bond is 



(4.1) 



in which Bq is given by (3.29) and 



IV^ 
lVs da ' 



da 



+ 



{T-tf 



B c + (T-t) 



93 
dr 



(4.2) 



We obtain {l\(z), IV£, IV2} from least-squares fitting, i.e. by minimizing 

n 

/X B ob B (^ S i) ~ ^modelC*) S H ^3 > lV 2)f \ 



(4.3) 



i=l 
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where B^ hs (t,Si) is the observed market price of a bond that matures at time Si and 
£^ odel (t, Si; IX, IV%, IV 2 ) is the corresponding model price obtained from (4.1). Here, n is 
the number of bonds that are traded at time t. For a fixed value of lX(z) it follows from (4.1 ) 
that {IV£,W$} can be determined as the least squares solution of 



9Bl(. q \ 1 

\ da \ L > p 



8BI ( Sl -t)* 
da "•" 2 



8B_l (S n -tf 



da 



B c + {S 1 -tf-§ 



B C + (S n -t)^ 



\ 



IVf 



(B^S^-B^SrJX)^ 
^B c ohs (t,S n )-B^t,S n ;fX)j 



Now, we vary l\{z) £ [0, Mi] and choose the point {IX, IV3, IV 2 f} that minimizes (4.3). Here, 
we take Mi = 1 guided by the results of Papageorgiou and Sircar (2008). 

Step 2. Estimation of {/, V{, V 2 £ , VI, Vg, V£, Vf} from the equity option data: These 
parameters are calibrated from the stock options data by a least-squares fit to the observed 
implied volatility. We choose the parameters to minimize 



i=l 



(Iobs(t, Ti, Ki) - Imodel(i) Ti, Kf, model parameters))' 



E 



(P obs (t, Tj, Ki) - P mode i(t, Ti, model parameters)) 2 
vega 2 (Ti,Ki) 



(4.4) 



in which I b s {t, Ti, Ki) and I mo de\(t, Ti, Kf, model parameters) are observed Black-Scholes im- 
plied volatility and model Black-Scholes implied volatility, respectively. The right hand side 
of (4.4) is from Cont and Tankov (2004), page 439. Here, P b S (t, Tj, ifj) is the market price 
of a European option (a put or a call) that matures at time Ti and with strike price Ki and 
-P mo dei(£> T, Ki; model parameters) is the corresponding model price which is obtained from 
(3.24). As in Cont and Tankov (2004), vega(Tj, Ki) is the market implied Black-Scholes vega. 



Let P (t,Ti,Ki;X(z)) be either of ( |3.44| and fl3.46[ ) with K = Ki and T = T t . Let us 
introduce the Greeks, 



91 



94 



-{T-t)x 2 

, d 3 Pn 



d 2 P 

dx 2 



dx 2 da ' 



95 



92 = —(T 
n 







2 d 2 P 



d 2 P n 



drjdx ' 
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1 



d 2 P n 



dadx da 



dPp | (T-t) 



dx I dx 2 



_ x d 2 Po 

dadx ' 
2 d 2 P 
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d_ 
da 



x- 



dPo 

dx 



Pc 



dx 2 



x 



97 = 

dPo 

dx 



(T-t) 2 2 d 2 P 



dx 2 



+ n 



dr 



dr 
(4.5) 



in which each term can be explicitly evaluated (see Appendix). 
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Now from (3.24) and the results of Section 



3.3 



(with I = 1), we can write 



^modelC*, Ti, Ki\ model parameters) = P Q (t, Tj, Kf, A) + V{gi(Ti, Kf, A) + V^giiTf Kf, A) 

+ VigaiTi, Kf, A) + Vtg^Ti, Kf, A) + Vig 5 (Ti, Kf, A) 
+ Vig^Ti, Kf, A) + Vtg 7 (T h Kf A) + V 2 5 g 8 (T u Kf, A). 

(4.6) 

First, let us fix the value of I. Then, from Step 1, we can infer the values of {A, V^ e ,V^}. 
Now the fitting problem in (4.4 ) is a linear least squares problem for { VJ, ^5 > ^6 > ^l*}- 



Next, we vary / G [0, 1] and choose {I, V{, V£, V£, V 6 e , Vf} so that (4.4) is minimized. 



4.3 Model implied CDS spread matches the observed CDS spread 

Let B c {t, T; I) denote the approximation for the price at time t of a defaultable bond that 



matures at time T, and has loss rate I (see (4.1 )). Let B(t, T) be the price of a risk-free bond. 



Then, the model implied CDS spead with maturity Tm is 



ds n t ^ B(t,T M ) - B c (t,T M ;l) , , 

Cmodell** J Mj = ^ ■ (4.7) 



Y,B c (t,T m ;l) 



m=l 



Recall that we have already estimated all of the model parameters in Section 4.2 using both 



corporate term structure data and the stock option implied volatility surface. Therefore, 



using (4.7) we can plot the model implied CDS spread over time and compare it with the 
CDS spread data available in the market. This is precisely what we do in Figures [T] and [2} We 
look at the time series c^ odel (t, 3) and c^ odel (i, 5) and compare them to the CDS spread time 
series of the Ford Motor Company. The match seems to be extremely good, which attests to 
the power of our modeling framework. 



By varying T]\j in (4.7) we can obtain the model implied term structure of the CDS 



spread. Figure [3] shows the range of shapes we can produce. 
4.4 Fitting Ford's implied volatility 

We will compare how well our model fits the implied volatility against the model of |Fouque 



et al. (2003), which does not account for the default risk and for the randomness of the 
interest rates. Although, we only calibrate seven parameters (hence we refer to our model 
as the 7-parameter model) to the option prices (see the second step of the estimation in 
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Section 4.2), we have many more parameters than the model of Fouque et al. (2003), which 



only has four parameters (we refer to this model as the 4-parameter model). Therefore, for 
a fair comparison, we also consider a model in which the volatility is a constant. In this 
case, as we shall see below, there are only three parameters to calibrate to the option prices, 
therefore we call it the 3-parameter model. 

Constant Volatility Model In this case, we take g\ = 02 = & m the expression for Pq 
in Corollary 3.1 . The expression for V6P . 1 remains the same as before. However, y/ePifi 
simplifies to 



ViPi,o = -(T 



x d 2 p | dP 

dadx da 



(4.8) 



This model has only three parameters, I, V{, Vf that need to be calibrated to the options 



prices, as opposed to the 4-parameter model of Fouque et al. (2003). 



As it can be seen from Figure [4] as expected our 7-parameter model outperforms the 
4-parameter model of Fouque et al.| ( |2003 ) as expected and fits the implied volatility data 
well. But, what is surprising is that the 3-parameter model, which does not account for the 
volatility but accounts for the default risk and stochastic interest rate, has almost the same 
performance as the 7-parameter model. 

The 7-parameter model has a very rich implied volatility surface structure, the surface has 



more curvature than that of the 4-parameter model of Fouque et al. (2003), whose volatility 



surface is more flat; see Figures [5] and [6] (The parameters to draw these figures are obtained 
by calibrating the models to the data implied volatility surface on June 8 2007.) The 7- 
parameter model has a recognizable skew even for longer maturities and has a much sharper 
skew for shorter maturities. 



4.5 Fitting the implied volatility of the index options 

The purpose of this section is to show the importance of accounting for stochastic interest 
rates in fitting the implied volatility surface. Interest rate changes should, indeed, be ac- 
counted for in pricing long maturity options. When we price index options, we set A = and 



our approximation in (3.24) simplifies to 

P € > 5 « P + yfeP lfl , (4.9) 
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in which Pq is given by Corollary 3.1 after settiing \(z) = 0, and 



1,0 



9 



Sx 2 2 dx 



dx 2 



, ire 2 d 3 P , T/£ 9 2 P ^ T/e ^0 



dx 2 da 



dr/dx 



dadx 
(4-10) 

Note that the difference of (4.9) with the model of Fouque et al. (2003) is that the latter 



allows for a slow evolving volatility factor to better match the implied volatility at the longer 



maturities. This was an improvement on the model of Fouque et al. (2000), which only has 



a fast scale component in the volatility model. We, on the other hand, by accounting for 
stochastic interest rates, capture the same performance by using only a fast scale volatility 
model. 



From Figure [7J we see that both (4.9) and Fouque et al. (2003) outperform the model of 
Fouque et al. ( |2000[ ), especially at the longer maturities (T = 9 months, 1 year, 1.5 years and 
2 years), and their performance is very similar. This observation emphasizes the importance 
of accounting for stochastic interest rates for long maturity contracts. 



Appendix: Explicit formulae for the Greeks in (4.5) 



When h{x) = {x — K) + , we can explicitly express the Greeks in (4.5) in terms of f(x 



h= exp(— x 2 /2) as 



2 d 2 C x/(di) 



x 



d ( 2 d 2 C 



dx 2 



dx \ dx 2 



1 



dx 



/V t ,T 



JL ( x ^ _ Co 

da \ dx 
d_ 
da 



-KB^T)[ T -^ + ^-^ T -^- 1 







2 



d 2 C 
dx 2 



(T - t exp(-/3(T -*))-! 



VtT 







2 



x _d_ (dCo^ = xfjdj) ( T-t + exp(-/?(r-t))-l 



dx \ da 
d ( dC 



dr V dx 



d (dC 



x 



+ 



dx \ dr) 
1 



<'i:r V 
C (l ) -KB c (t,T) 

1 If] 



l-exp(-/3(T-t)) 



N{d 2 



N{d 2 



f(d 2 ) 

JVtT 



xf(di 



^ v ^ (T "' )+ ^ (exp( " /?(T " t)) " 1 
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3 (exp(-2/3(T-t))-l)) 



log 



-3/2 



KB? T *> T 



1 



^■JvCf 



'2pi<j 2 2rj\ f2 Pl a 2 4rA rj (2 Pl a 2 3r] 

-y~ + p2 ) ( T ~ *) + ( —pT + J e M-P(. T ~ *)) - exp(-2^(T - t)) - \ — ^ + 



2 



13* 
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Figure 1: Ford 3 year CDS annual spread time series from 9/18/2006-3/13/2007. 

Spread implied by model is pink solid line, real quoted spread is blue broken line. Ford's 3 

year CDS spread time series is not available from 3/13/2007 until 8/31/07 in our data source. 
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Figure 2: Ford 5 year CDS annual spread time series from 9/18/2006-6/8/2007. 
Spread implied by model is pink solid line, real quoted spread is blue broken line. 
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Figure 3: CDS Term Structures (4.7) can produce: 
Legend 

-x-, blue (The parameters are obtained from calibration to 11/13/2006): a=0.0037, /3=0.0872 
r] = 0.0001, r = 0.0516, /(loss rate) = 0.283, X(z) = 0.0459, [V£,Y$] = [0.0425,0.0036]. 

-squares-, black (The parameters correspond to 6/18/2006): a = 0.0045, [3 = 0.0983, 
7] = 0.0002, r = 0.0516, I = 1, A = 0.012, [V£, Y$] = [0.0185,0.0025], 

-diamonds-, red (The parameters correspond to 9/22/2006): a = 0.0039, (3 = 0.0817, 
7] = 0.0012, r = 0.0496, 1=1, A» = 0.017, [V£,Vf] = [0.0067,0.0005] 
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Figure 4: Implied volatility fit to the Ford call option data with maturities of 
[17,45,72,168,285,643] calender days on April 4, 2007. 

Model is calibrated aross all maturities but we plotted the implied volatilities for each matu- 
rity, separately. Here, stock price (x) = 8.04, historical volatility (02) = 0.3827, one month 
treasury rate (r) = 0.0516, estimated correlation between risk-free spot rate(one month trea- 
sury) and stock price (pi) = -0.0327. Also a = 0.0037, (3 = 0.0872, r) = 0.0001 which are 
obtained with a least-square fitting to the Treasury yield curve on the 4th of April. 
Legend: 

'o', empty circles = observed data; 

'x', green = stochastic vol+stochastic hazard rate+stochastic interest rate = the 7-parameter 
model; 

small full circle, blue = constant vol+stochastic hazard rate+ stochastic interest rate = the 
3-parameter model 

red = The model of Fouque et al. (2003) which has constant interest rate+stochastic vol 
(slow and fast scales) = the 4 parameter model. 
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Figure 5: Implied volatility surface corresponding to (4.6), the 7-parameter model. 
Here, a = 0.0063 1 p = 0.1034, rj = 0.012, r 
a 2 0.2576, \(z) 0.027, (Vf, V£, Vf, V£, V£, V{, V 2 S ) 

(0.9960, -0.0014, 0.0009, 0.0104, -0.6514, 0.3340, -0.1837, -0.0001). 
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Figure 6: Implied Volatility Surface corresponding to the 4-parameter model of Fouque et al. 
(20031. 



Here, 



0.046, 



average volatility=0.2546, 
(2003[) are choosen to be (V£, Vg, V$ , 



and the parameters in (4.3) of Fouque et al. 
0.0164 



:! . .„..;, - (-0.0164,-0.1718,0.0006,0.0630). Note that 
the parameters here and Figure [5] are both obtained by calibrating the models to the data 
implied volatility surface of Ford Motor Company on June 8, 2007. 
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Figure 7: The fit to the Implied Volatility Surface of SPX on June 8, 2007 with maturities 



[30,60,91,122,152,182,273,365,547,730] calender days. Recall from Section [471] that we use 
standardized options from the OptionMetrics. 

Models are calibrated aross all maturities, but we plot the implied volatility fits separately. 
The parameters are: stock price (x) = 1507.67, dividend rate = 0.0190422, historical volatility 
((T2) = 0.1124, one month treasury rate (r) = 0.0476, estimated correlation between risk-free 
spot rate(one month treasury) and stock price (pi) = 0.020454. Also, a = 0.0078, = 0.1173, 
r) = 0.0241, which are obtained from a least-square fitting to the Treasury yield curve of the 
same day. 
Legend 

'o', empty cirles = observed data, 



'x", green = Implied volatility of (4.9), 

'*', red = Implied volatility oflFouque et al.|(|2003|) 



small full circle, blue = Implied volatility of Fouque et al. (2000). 
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